Phys. Rev. D 86, 044021 (2012) 



AEI-2012-034 



larXiv: 1204 .25501 



Diffusion in quantum geometry 



(N 

O 

(N 

< 



Oh- 

CD ■ 

> 
O 

in 
o 



X 



Gianluca Calcagni 
Max Planck Institute for Gravitational Physics (Albert Einstein Institute), 

(Dated: April 11, 2012) 



Am Muhlenberg 1, D-14-i76 Golm, Germany 



The change of the effective dimension of spacetime with the probed scale is a universal phe- 
nomenon shared by independent models of quantum gravity. Using tools of probability theory and 
multifractal geometry, we show how dimensional flow is controlled by a multiscale fractional diffu- 
sion equation, and physically interpreted as a composite stochastic process. The simplest example 
is a fractional telegraph process, describing quantum spacetimes with a spectral dimension equal to 
2 in the ultraviolet and monotonically rising to 4 towards the infrared. The general profile of the 
spectral dimension of the recently introduced multifractional spaces is constructed for the first time. 
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The spectral properties of effective quantum geome- 
tries show that the ultraviolet (UV) finiteness of indepen- 
dent theories of quantum gravity is universally associated 
with a lower spectral dimension of spacetime (typically, 
ds ~ 2) at small scales, while dg ~ 4 in the infrared (IR). 
Instances are causal dynamical triangulations (CDT) y|, 
asymptotic safety (QEG) 0, 3, spin foams 3,01: non- 
commutative geometry jg] , Hofava-Lifshitz gravity 0] , 
and other approaches [8| . The change of dimension with 
the probed scale is known as dimensional reduction or di- 
mensional flow Understanding its physical meaning 
is an important piece of the puzzle of quantum grav- 
ity, since multiscale behavior is deeply related to the 
renormalization properties of these theories. Differen- 
tial geometry and ordinary calculus, as employed in gen- 
eral relativity and field theory, are inadequate to study 
this and other properties of quantum spacetimes, and 
stochastic processes and multifractal geometry can offer 
powerful tools of analysis and novel insight. While there 
is the tendency to label all multiscale spaces as "frac- 
tal," the accumulated knowledge from these branches of 
mathematics permit to make sharper statements about 
the geometric and physical properties of quantum-gravity 
models. This philosophy inspired the revisiting of a re- 
cent problem, the construction of quantum field theories 
in fractal spacetimes, under a fresh perspective focused 
on an effective continuum geometry [lO| , in particular via 



the formalism of multifractional spacetimes [11 



After a sketch of the classical situation, we will argue 
that quantum geometry effectively modifies the diffusion 
equation. A critical appraisal of the latter will allow us to 
classify quantum geometries in terms of stochastic pro- 
cesses on one hand, and to get a precise back-up to the 
notion of "fractal spacetime" on the other hand. The 
aim is to reexamine the spectral dimension starting from 
its foundation and provide a general, model-independent 
and analytic description of dimensional flow, conflrmed 
by quantum-gravity examples. This is possible thanks to 
the presence of universal features in the flow [T^] . 



P of the transport equation must be nonnegative at all 
points. If P is normalized to 1, it is interpreted as the 
probability to find the diffusing test particle at a given 
point. This probability distribution describes a stochas- 
tic process, i.e., a sequence or collection of random vari- 
ables. We shall use the condition P > as one of the 
guiding principles to identify the random process associ- 
ated with a given behavior of quantum geometry. Here 
we do not pay attention to the techniques solving the 
diffusion equations; an expanded discussion is in [12|. 

Classical spacetimes. In a smooth classical spacetime 
with D topological dimensions, the diffusion equation is 



{d„~yl)P{x,x\a) = Q. 



(1) 



The parameter cr > acts as an abstract "time" vari- 
able via the diffusion operator d„ , an ordinary first-order 
derivative. Writing a = I'^a in terms of a length scale I 
and a dimcnsionless parameter ct, Eq. ([T|) is recast in the 
form {da — PVl.)P = 0. The spatial generator is the 
Laplacian in the given metric background in Euclidean 
signature, x' is the initial point where diffusion starts. 
In translation-invariant spacetimes, the heat kernel P de- 
pends on the difference x — x' , but in fractional spaces 
with nontrivial measure this is no longer true; therefore 
we keep the notation P{x, x' , a) separate from the often- 
employed u{x,a) (fixing x' = 0). Equation ([T|) is not 
completely specified without the set of initial conditions 
at (T = 0. The choice P(x,x',0) = 5{x — x') describes 
diffusion of a point particle starting at x = x' . Extended 
shapes of the probe are possible [8|, but the pointwise 
one allows us to explore the local structure of spacetime. 

Ignoring curvature, the solution P must be non- 
negative for all X and x' , and normalized as 
J d^xP{x,x' — 1. From the spatial trace 
of P, one gets the return probability P((t) := 
(/ dfx)^^j d^xP{x,x,a) and the spectral dimension 



^dlnP(cr) 
diner 



(2) 



For a diffusion process to be meaningful, the solution When divergent, the volume prefactor in the definition 



2 



of V can be regularized. In the case of a translation- 
invariant background, it cancels out with the position 
dependence in the numerator and 'P{a) = u{Q,(j). 

The normalized solution of Eq. ^ is the Gaussian 
heat kernel P(x,x',ct) = wi(r,cr) e-""^ 1'^^"^ j (\T:a)^ 1"^ , 
where Clearly, P > 0. The return 

probability and spectral dimension read V oc cr^^/^ and 
ds = I?, respectively. There is no quantitative distinction 
between spectral and topological dimension. They are 
also equal to the Hausdorff dimension dn of spacetime, 
determining the scaling law of the volume of a D-ball of 
radius i?, V*-^-* oc R^^ . Notice that, because I is the only 
scale, it is not possible to define a hierarchy of scales and 
the geometry (and ds) is scale independent. 

The ordinary diffusion equation ([T]) is associated with 
a Wiener process B{a), the standard Brownian motion. 
B is continuous in a with probability 1, B{0) = x' , and 
the increments of B are governed by the Gaussian distri- 
bution ui, B{a) — B{a') ^ wi(0,(T — a') for a' < a. 

Quantum geometry with fixed dimension. Quantum ge- 
ometry can emerge either by definition of a nonstandard 
texture of spacetime [6l-[8l [Tcil . [T]| or from the quantization 
of gravity or for both reasons. Since the spectral di- 
mension ds becomes then anomalous, quantum geometry 
effectively modifies Eq. ([T]). Modifications affect either 
the initial condition P(a;, a;',0) (by a change of the defi- 
nition of "point particle" on a quantum manifold or of 
the delta distribution as in multifractional spaces [Tl|). 
or the operator dc as in QEG or the Laplacian 
(via the change in the differential structure and/or the 
presence of one or more fundamental quantum scales £„, 
such as the Planck scale or the label-dependent lengths 
of the simplices in a cellular complex), as in QEG Hi, 
Hofava gravity 0, multifractional theory [fli. and CDT, 
spin foams and simplicial gravity in general [ll, |3j ■ 

Dimensional flow is still inadequately understood and 
a classification of the possible diffusion equations should 
help to control the physics of the above-mentioned (as 
well as other) models of quantum gravity. It is instruc- 
tive to specialize first to the case of fixed dimensional- 
ity (no scale hierarchy). We concentrate on the contin- 
uum formulation of fractional calculus, which guarantees 
anomalous (in particular, fractal) geometric properties 
of spacetimes [Tl|i Ell ^"^^ anomalous correlations in dif- 
fusion problems (e.g., (l^). Multifractional theory is a 
model of quantum gravity in its own right (like Hofava 
gravity, it is a traditional perturbative field theory but 
built on an anomalous spacetime), although it can serve 
as a framework to understand other proposals 0, Q . We 
ignore curvature. The latter modifies the spectral prop- 
erties of spacetime even in a classical setting, except in 
the UV limit ct — > 0. Quantum effects, however, often 
modify spacetime globally even in the absence of curva- 
ture, which motivates the assumption (see also [1, Q). 
We replace 9^ with the operator 

1 Hy,{oodf^ + oodf^) r , — ^1 

lCl,a-= -===—^ r \JVo,[x)- , 3) 

^j%Ax) 2cos(7r7) L J 



where Va{x) oc ja^/^I" ^ is the measure weight of 
the ambient space (the singularity in a;^ = is in- 
tegrable and does not pose particular problems for 
the classical and quantum dynamics), < a < 1 
and < 7 < m are real parameters, and we make 
use of left and right Liouville-Caputo fractional deriva- 
tives: for each direction, {ood"^^ f)ix) oc J^^dx'{x — 

x'r-^-^w^f{x'), Udy){x) oc {~iy" jfdx'{x' ~ 

x)^-^-^i_d^f[x'). When 27 = m is integer, oo^" = 
(-l)"ooi9"' = (9". Definition © is such that, in a suit- 
able domain, the operator /C^.q is self-adjoint jTll [l^. 

We classify the stochastic and geometric properties as- 
sociated with the diffusion equation 

(9. - /C^,a) -P = , (4) 

with initial condition P I o-=o = \vn{x)Vn {x')]~^^'^5{x~x'). 
One can show that ds = Da/^ |12| . 

(i) When 7 = 1 = 0;, we recover ordinary diffusion 
and ds = D. For a 7^ 1, this is ordinary Brownian 
motion but on a fractal spacetime with ds = Da- 

(ii) For < 7 < 1 and a — I, we have a supcrdiffusive 
Levy process with ds = D / j > D = da- This does 
not correspond to a fractal spacetime (ds < dn for 
fractals). For < 7, a < 1, one has a Levy process 
on an anomalous spacetime (fractal if a < 7). 

(iii) When 7 > 1, the solution of Q is no longer nonneg- 
ative definite and the equation must be modified. 
In fact, one can include a source term S{x,x',a), 
which does not alter the spectral dimension. Hence, 
overlooking the check that P > for the Ansatz (|4]) 
might result in the correct spectral dimension but 
a wrong diffusion equation. 

Processes associated with nonhomogeneous equations 
may be non-Markovian (i.e., future states depend also 
on past states) even if they are meaningful in a proba- 
bilistic sense. An example is the quartic equation 

(a, - V^) u{x, a) = {nar^/^\/lu{x, 0) , (5) 

with source given by the initial condition. The solu- 
tion gives the same ds = D/2 a.s the naive Eq. with 
7 = 2, but the presence of the source guarantees that 
M > 0. Equation ([5]) governs an iterated Brownian mo- 
tion (IBM) Given two independent Wiener pro- 
cesses Pi, 2, IBM is defined as Xibm(o') Bi[\B2{<t)\], 
where B2 acts as a clock to Pi. In general, there exists a 
deep connection between higher-order diffusion equations 
with integer time, iterated stochastic processes, and dif- 
fusion equations with fractional time (5^ — V^)w = 0, 
with < /3 < 1 and the left Caputo derivative (lower 
terminal ct' = — 00 replaced by a' = 0). The solution u is 
positive definite • The process described by this equa- 
tion (fractional Brownian motion) is subdiffusive and on 
average it takes longer (with respect to Brownian motion) 
for the particle to cover a certain distance. Equation ([5]) 
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can be regarded as the "iteration" of the fractional equa- 
tion with j3 = 1/2, with the same solution u. 

Extending the discussion to a nontrivial spacetime 
measure, the spectral dimension associated with the frac- 
tional diffusion equation {d^ — /C^^q,)P = is ds = 
(/3/7)(iH, where dn = Da [HI A fractal configu- 

ration is obtained whenever /? < 7. 

Multiscale quantum spacetimes. We now make a con- 
ceptual step of relevance for quantum gravity. In complex 
systems, the adoption of a diffusion equation is motivated 
by phenomenology. Given a set of experiments evidenc- 
ing some anomalous scaling laws, one proposes an ad-hoc 
diffusion equation reproducing those scalings. The the- 
ory is then further tested against experiments. Or else, 
one defines the stochastic process underlying a certain 
physical system, and from its probability distribution one 
infers the correct diffusion equation. For instance, IBM 
provides a stochastic description of diffusion of a particle 
trapped in a crack [l6| , the latter being modeled by a ran- 
dom fractal set whose pattern is the graph of a Brownian 
motion. In quantum gravity, on the other hand, we do 
not have experiments but fragmentary knowledge such as 
the existence of anomalous scaling behaviors in the UV. 
This information determines the differential order of the 
operators in the diffusion equation, i.e., the number and 
values of the plateaux in the dg profile, but it may be 
unable to fix the diffusion equation univocally. In par- 
ticular, the monotonic slopes between the plateaux can 
vary from model to model of the same physical system, 
but they may be not falsifiablc features; their nonunique- 
ness can be usually traced back to details of the theory 
such as regularization schemes [1, [l^]- Also, dimensional 
flow may be insensitive of the presence of source terms 
[Eq. ^ without source or with flipped sign in front of 
would still give the same dg], and we must resort to 
positivity of the solution P to fix more details of the dif- 
fusion equation. In turn, once we determine a reasonable 
diffusion equation with probabilistic interpretation, we 
can also find the stochastic process associated with that, 
thus physically characterizing quantum geometry. 

Without further input from the theory except the 
UV and IR behaviors, we can reproduce the whole di- 
mensional flow by applying the techniques of multiscale 
phenomena and multifractal geometry to the texture of 
spacetime itself [l^l ■ The generalization of the diffusion 
equation (with solution P > 0) to a multiscale process is 
realized by summing over all possible values of a, {3, 7: 



(?n9^" " Cn/C^„,a„) P{x, x' , a) = S{x, x\ Cr) , (6) 



where ^„ and C„ arc dimcnsionful couplings which de- 
pend on the characteristic scales of the system. Typi- 
cally, there is only a finite number N of terms in phys- 
ical systems, so the sum representation © is realistic. 
The number TV — 1 of characteristic scales (hidden in ^ 
and C) determines the number TV of plateaux (asymptotic 
regimes) in the profile of c?s- A multiscale phenomenon 
is always defined by the relative size of the scales, not by 
an absolute hierarchy. This means that we can choose 



any of the TV scales £„ to represent the scale £ probed 
by a measurement. If we order the scales of the sys- 
tem as £i < I2 < ■ ■ ■ < ^AT, we can take the largest as 
£ = £n- Thus, there are TV — 1 (not TV) scales with the 
physical meaning of characteristic lengths. The spectral 
dimension is fixed when TV = 1; for TV = 2, it has two 
asymptotic values dg ^ c?si.2 in the regimes £ <^ £1 and 
£ £1, with a monotonic transient phase in between; for 
TV = 3, there will be also an intermediate plateau where 
ds dg3; and so on. 

The first example is an interaction of Gaussian and 
anomalous dynamics which can describe certain turbu- 
lent media The diffusion equation is (da ~ — 
CiA^7,i)m = 0, u(x,0) = 5{x), where < 7 < 1 and 

we write the constant Ci = ^1 ''^ in terms of a char- 
acteristic length. As the analytic solution shows, the 
transport is of Levy type at large scales £ = ^ £1 
(ds ^ 1/7 > 1 = dn) and normal at small scales £ ^ £1 
(ds '--^ 1). From the perspective of quantum spacetimes, 
this model is multiscale but not multifractal. Profiles 
of ds overshooting the Hausdorff and topological dimen- 



sions appear also in lattice 



and noncommutative ge- 



ometries (last reference in [6|), with some caveats [l2l |. 

A second example is a fractional diffusion equation 
with two diffusion operators da^'^ . To see its neat 
stochastic interpretation, we recall some results on the 
so-called telegraph processes [l^. A telegraph process 
is defined as V{a) = V{Q) (-l)^('"), where V{a) is the 
velocity of a particle at time a running on the real line, 
V{0) is the initial velocity which is ±c with equal prob- 
ability, and M is the cumulative number of events of 
a homogeneous Poisson process with rate A > 0. The 
velocity of the particle flips direction back and forth, 
hence the name "telegraph." The position of the parti- 
cle at time a is the integrated telegraph process T{a) ~ 
V{0) ds We consider a composite process 

called l3rownian-time telegraph process or fractional tele- 
graph process^ Xftp(o') := T[|i3(cr)|]. This motion is gov- 
erned by the diffusion equation [d„+2\dy^ — c^&l)u = 0. 
The solution of the fractional telegraph equation and its 
generalizations is nonnegative and unique [l^. In the 
double limit A,c — !• 4-00, \/(? — const, the stochas- 
tic process reduces to an IBM. Recasting these results 
in the language of multifractal spacetimes and extending 
to D dimensions, we set [a] = 0, c = £j as the charac- 
teristic scale, and £ ~ £*/(2A) as the probed scale. In 
the limit £ 3> diffusion is Gaussian (Brownian pro- 
cess, ds ^ D). At small scales i <^ £^:, one reaches a 
regime where diffusion is fractional and given by an IBM 
(ds D/2). In between, diffusion in quantum spacetime 
obeys the law of a fractional telegraph process. 

The monotonic profile ds{£) of this single-scale space- 
time can be plotted from the analytic form of the return 
probability [12| . The probability distribution for more 
complicated multiscale spacetimes can be computed as 
well, but here we show how all these proflles are easily 
reproduced in the framework of multifractional geometry 
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FIG. 1. The spectral dimension ds{£) in D = 4 for a multi- 
fractional model and normal diffusion (/? = 1 = 7, dg = dn) 
with single scale (dashed curve) and two scales (solid curve). 



(q!„ 7^ 1) when /S = 
coefficients <^„ in ([ 



1=7 and S ~ 0. Wc argue that the 
|) have the natural form Cw = 1 and 



(7) 



where t = In ■ First, we notice that the Laplacians all 
have the same order 2, so the coefficients C„ all have the 
same scaling dimension and, in particular, we can always 
make them dimensionless. Write C„ as the ratio of some 
length scales, Cn = {I A,n / 1 b ,n)'^ ■ Without loss of general- 
ity, one can choose q = 2 so that the spatial generator of 
the diffusion equation can be rendered dimensionless, in 
the form X]n(^^,")^^i,a„ ■ Now, the nth term dominates 
over the others at scales i <^ £n, while at scales smaller 
than £n-i the (n— l)th term takes the lead, so the small- 
est possible scale i at which the nth term dominates is 

£ ^ £n-l- Thus, we set lA,n = £n and lB,n = £ — £n-l- 

In other words, the dimensional flow is always measured 



starting from the lowest of two scales £n-i to the next 
£n, and relatively to the latter, which sets a gauge for the 
rods. Since £ = £n is the probed scale. Cat = 1- 

We can plot the spectral dimension for any given profile 
a{£). Upgrading on plj . we motivate a realistic profile 
a{£) as an approximation of the sum in ([6]). Consider 
first the N = 2 case with ai ^ 1 and 02 = 1. In one 
dimension, and by Eq. ©, {d^ + CiK^f a JP = {(1 + 
Ci)/CL^(,)+Ci(l-«i)V[4(l + Ci)a^']m where ai{£) := 

[1 + Ci(^)ai]/[1 + Ci(^)]. For both smaU and large Ci 
the kinetic term in this expression dominates over the 
potential term, so the profile ai{£) defines an effective 
fractional charge throughout the dimensional flow. With 
N coefficients a„, = 1, the effective charge reads 



1 



aN-i{i) ■■= 



1 + EL~i'C«W 



(8) 



For two entries {N = 2, ai = 2/D, a2 = 1, one 
scale), dimensional flow is such that ds ^ D in the IR 
and dg ^ Dai = 2 in the UV, with no intermediate 
regime in between. This is the type of flow considered in 
[lol . [Til and is shown in Fig. [1] (dashed curve), in quan- 
titative agreement with the fractional-telegraph profile 
which is matched by tuning the scale length [3[. A two- 
scale profile ds{£) = ^a2{£) with ai = 1/2, a2 = 1/3 
and £2 = 10£i is also plotted (solid curve). At ^ = 0, 
ds = 2. At £ ^ £1, the spectral dimension acquires the 
minimum value = 4/3. At scales £ < £2, the diffusion 
process corresponds to a recurrent random walk, where 
dg < 2. Well above the larger critical scale, £ ^ £2, 
ds hits the IR value ~ 4. Notably, the profiles in Fig. 
[1] reproduce the dimensional flow of asymptotically-safe 
quantum gravity in two different realizations, without or 
with matter. This is consistent both with the fact that 
the single-scale curve of QEG does come from a fractional 
telegraph process equation and with a reinterpretation of 
the renormalization group flow in terms of measurements 
in multifractional geometry Q . 



[1] J. Ambj0rn, J. Jurkiewicz, and R. Loll, |Phys. Rev.| 
Lett. 95, 171301 (2005); D. Benedetti and J. Henson, 
|Phys. Rev. D 80, 124036 (2009), 
[2] O. Lauscher arid M. Renter, |J. High Energy Phys.| 
[10 (2005) 050; M. Reuter and F. Saueressig, 



J. High Energ y Phys. 12 (2011) 012 
[3] G. Calcagni, A. Eichhorn, and F. Saueressig (work in 
progress) . 

[4] L. Modesto, [Classical Quantum Gravity 26, 242002] 
•^09); F. CaraveUi and L. Modesto. [arXiv: 0905 .21701 
E. Magliaro, C. Perini, and L. Modesto, 
arXiv: 091 1701371 

[5] G. Calcagni, D. Oriti, and J. Thiirigen, 
larXiv: 1208 .03541 G. Calcagni, D. Oriti, and J. 
Thiirigen (wo rk in progress) . 

[6] A. Connes, [J. High Energy Phys. 11 (2006) 081[ D. 



Benedetti, [Phys. Rev. Lett. 102, 111303 (2009)[ M. 
Arzano, G. Calcagni, D! Oriti, and M. Scalisi, 
Phys. Rev. D 84, 125002 (2011)[ E. Alesci and M. 



Arzano, [Phys. Let t. B 707, 2 72 (2012) [ 

[7] P. Hofava, |Phys. Rev. Lett. 102, 1613 01 (2009)'; 
T.P. Sotiriou, M. Visser, and S. Weinfurtner, 
[Phys. Rev. Lett. 107, 1313 03 (2011) 
[8] L. Modesto and P. Nicolini, (Phys. Rev. D 81, [ 



104040 (2010) E. Spallucci, A. Smailagic, and P. 



Nicolini, [Phys. Rev. D||73, 084004 (2006)[ L. Modesto, 



larXiv: 1107.24031 

[9] G. 't Hooft, in Salamfestschrift, edited by A. 
All, J. Ellis, and S. Randjbar-Daemi (World Scien- 
tific, Singapore, 1993); S. Carlip, larXiv: 0909 .33291 
[iEXiv : 1009JJ^ T.P. Sotiriou, M. Visser , and S. We- 
infurtner, [Phys. Rev. D 84, 104018 (2011)[ 



5 



[10] a Calcagni, Phys. Rev. Lett. 104, 251301 (2010) 



J. High Energy Phys. 03 (2010) 120.: Phys. Lett. B 697,, 



251 (201i: 



[11] G. Calcagni, 



Phys. Rev. D 84, 061501(R) (2011) 
I arXiv: 1106.57871 [Adv. Theor. Math. Phys.] ~(to 
be pubhshed); J. High Energy Phys. 01 (2012) 065 



G. Calcagni and G. Nardelh, arXiv : 1 202 . 53831 

[12] G. Calcagni, arXiv : 1205 . 504Q 

[13] G.M. Zaslavsky, Phys. Rep. 371, 461 (2002)', A.A. Kil- 
bas, H.M. Srivastava, and J.J. Trujillo, Theory and Ap- 
plications of Fractional Differential Equations (Elsevier, 
Amsterdam, 2006). 
[14] H. AUouba and W. Zheng, |Ann. Probab. 29, 1780| 



(2001); R.D. DeBlassie, "Ann. Appl. Proba b. 14, 1529| 
(2004); B. Baeumer, M.M. Meerschaert, and E . 
Nane, [Trans. Am. Math. Soc. 361, 3915 (2009) 



Beghin, E. Orsingher, and L. Sakhno, 



[15] 



]Stoch. Anal. Appl. 29, 551 (2011) 




E. Orsingher and L. Beghin, 


Ann. Probab. 37, 206| 



(2009) 



[16] K. Burdzy and D. Khoshnevisan, [Ann. Appl. Probab. 8,| 
|708 (1998)] 



[17] H. Weitzner and G.M. Zaslavsky, [Chaos 11, 384 (2001) 
[18] K.J. Hochberg and E. Orsingher, 'J. Theor. Probab 
19, 511 (19 96T1 E. Orsingher and L. Beghin, 
|ProEabrTheory|Relat. Fields 128, 141 (2004) | 



